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Stability Analysis of Integral Delay Systems with Mnltiple Delays 


Bin Zhou* Zhao-Yan Li^ 


Abstract 

This note is concerned with stability analysis of integral delay systems with multiple delays. To 
study this problem, the well-known Jensen inequality is generalized to the case of multiple terms by 
introducing an individual slack weighting matrix for each term, which can be optimized to reduce the 
conservatism. With the help of the multiple Jensen inequalities and by developing a novel linearizing 
technique, two novel Lyapunov functional based approaches are established to obtain sufficient stability 
conditions expressed by linear matrix inequalities (LMIs). It is shown that these new conditions are 
always less conservative than the existing ones. Moreover, by the positive operator theory, a single LMI 
based condition and a spectral radius based condition are obtained based on an existing sufficient stability 
condition expressed by coupled LMIs. A numerical example illustrates the effectiveness of the proposed 
approaches. 

Keywords: Stability of integral delay systems, Multiple Jensen inequality, Linearization technique. 
Positive operator theory, Spectral radius 


1 Introduction 

An integral delay system (IDS) in the form of x (t) = /A F (s) a; (t -|- s) ds, where F (s) is a matrix function 
with bounded variation, has many important applications in the study of time-delay systems (see, for ex¬ 
ample, [3], 0120] and ED). In Hale’s book 0 this class of IDSs were named as D operators, which were 
also treated as generalized difference equations there, and their stability is necessary for the stability of the 
associated neutral time-delay systems. This class of IDSs also come from the model reduction approach for 
stability analysis of time-delay systems, which are frequently named as the additional dynamics (see 0 and 
[ini) and their stability is necessary for the stability of the transformed time-delay systems. This class of 
IDSs are very closely related with the predictor feedback control of linear systems with input delays, for 
example, 


• In 0 , the author proved that the numerical implementation of the predictor feedback for linear systems 
with input delays is safe only if an IDS is exponentially stable. 

• In [14] it is shown that the stability of this class of IDSs is necessary for the robust stability of linear 
systems with input delay by the well-known predictor feedback. 

• In [12] we have shown that an input delayed linear system by the so-called pseudo-predictor feedback 
is exponentially stable if and only if an IDS is exponentially stable. 


Stability analysis of IDSs can be traced back at least to Cruz and Hale 0 , Henry 0 and Melvin m, in the 
study of the stability of neutral time-delay systems 0. When the right hand side of the IDS only contains 
terms at some time points, a general theory was build in 0. This class of IDSs have received renewed interest 
in recent years. A general stability theorem was build in m and was later applied on different forms of 
IDSs (see [T3], [S], [IH], and the references therein). In the paper 0 stability conditions are derived for IDS 
with matrix discrete-continuous measures. 
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^Zhao-Yan Li is with the Department of Mathematics, Harbin Institute of Technology, Harbin, 150001, China. Email: 
zhaoyanlee@gmail.com; lizhaoyan@hit.edu.cn. 
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In this note we also study stability analysis of this class of IDSs in which F (s) is a piecewise constant 
matrix function, namely, the system can be expressed as x {t) = I-t 2 ; (t + s) ds where Ai and 

Ti are constants (r^ can be unknown, time-varying, yet bounded). This class of IDSs with multiple time 
delays have been investigated in m with the help of the well-known discrete-time and continuous-time 
Jensen inequalities. In this note, by recognizing that the jointed usage of the conventional discrete-time and 
continuous-time Jensen inequalities requires that all the integration terms (t + s) AjQAiX {t -|- s) ds 

share the same weighting matrix Q, we first establish a so-called multiple Jensen inequality, by which as 
well as some novel Lyapunov functionals and a new linearization technique, every individual integration 
term possesses a different weighting matrix, which can introduce more weighting matrices that will be 
optimized to reduce the conservatism of the corresponding sufficient stability conditions. Indeed, it is shown 
in both theory and by numerical examples that the stability conditions obtained by the multiple Jensen 
inequalities are always less conservative than that obtained by the jointed usage of the conventional discrete¬ 
time and continuous-time Jensen inequalities. Another contribution of this note is that, with the help of 
the positive operator theory, we are able to establish an equivalent linear matrix inequality (LMI) based 
stability condition involving only one constraint and one decision matrix and an equivalent spectral radius 
based stability condition of some existing stability conditions that are expressed by a set of coupled LMIs. 
These two equivalent stability conditions are appealing in both theory and in computation. 

The remaining of this note is organized as follows. The problem formulation and some preliminaries are 
given in Section [2] In Section |3l two kinds of LMIs based sufficient conditions are established with the help 
of the multiple Jensen inequalities, some novel Lyapunov functionals and a new linearization technique. A 
spectral radius based stability condition and a single LMI based stability condition are then established in 
Section |4] in which a comparison of the proposed approaches and the existing one will also be carried out. A 
numerical example is worked out in Section [5] to illustrate the effectiveness of the proposed approaches and 
finally Section |6] concludes the note. 


2 Problem Formulation and Preliminaries 


Consider the following integral delay system (IDS) with multiple delays 


N 

x{t) = Y^ A^ 


X (t -h s) ds, 


( 1 ) 


where Ai G R”^", i G I[I,A^] = {1,2,..., A^}, are given matrices and Ti,i G I [1, TV], are given scalars and 
are such that 

0 < Ti < T = max {xi}, T € I [I, TV]. (2) 

zei[i,iv] 

Let ip G ^n,T be an initial condition for m and x{t) = X (T, p),'it > 0 be the corresponding solu¬ 
tion of o satisfying x{t) = p{t),'it G [—t, 0). Here '^n,T denotes the Banach space of continuous vec¬ 
tor functions mapping the interval [—t, 0] into R” with the topology of uniform convergence. We say 
that the IDS o is exponentially stable if there exist two positive constants a and /3 such that ||a;(T)|| < 
asup^g[_^^o) llT5(s)||e-'^*,Vt > 0. 

The IDS ([II arises when some transformations are made on differential-difference systems [1^. In this note 
we are concerned with the stability analysis of the IDS ([T]). By choosing some suitable Lyapunov functionals 
and developing a new linearization technique for handling nonlinear matrix inequalities, we will establish two 
classes of LMIs based sufficient conditions guaranteeing the exponential stability of the IDS ([1]). Moreover, 
with the help of the positive operator theory, we will also provide a spectral radius based sufficient stability 
condition. The relationships among these different sufficient conditions are also revealed. Our results improve 
those in [16) . Both theoretical analysis and numerical examples will demonstrate that the obtained results 
are always less conservative and more efficient than the existing ones especially those in [16] . 

The following general Lyapunov stability theorem for the IDS ([T|) will be used later in this note. 


Lemma 1 Udj/ The IDS m 


is exponentially stable if there exists a differentiable functional V : 'ffn,T R 


2 


( 3 ) 


and three positive constants on,i = 1 , 2 , 3 , such that 


/ O pO 

\\x {t + 9)\\^ dO < V {xt) < a 2 J \\x{t + 9)\\^A9, 


Vixt)<-a 3 \\x{t + 9)\\'^ d9. 


(4) 


At the end of this section, we give the following technical lemma which is helpful for the linearization of 
nonlinear matrix inequalities in the sequel. 


Lemma 2 Let S G R.'ix" Q g he two positive definite matrices. Then Q < S ^ if and only if 

there exists a matrix R G R^xn such that 

R^QR +S - {R +R^) <0. (5) 

The same statements hold true if “< ” in the above two inequalities are replaced by “< ”. 


Proof. It follows from ([S]) that R is nonsingular. Since {R — S)^ S ^ {R — S) >Q, namely, 

-R^S-^R <S-{R + R^), 

we have from ([S]) that 

Q < -R-'^iS -{R + R^))R-^<R-'^R'^S-^RR-^ = R-fi 
On the other hand, if Q < S~^ is satisfied, then ([S]) is satisfied by choosing R = S. ■ 


( 6 ) 

(7) 


3 The Multiple Jensen Inequality Based Stability Conditions 

3.1 The Multiple Jensen Ineqnality 

We first recall the following well-known Jensen inequality. 


Lemma 3 For any positive definite matrix Q > 0, a positive number r > 0, and a vector valued function 
oj : [—T, 0] —>■ R” such that the integrals in the following are well-defined, then 


J uj{s)ds'^ ^ (y* uj{s)ds^ • 

Moreover, for a series of vectors z € I [1, A^], there holds 


N 


N 


N 




( 8 ) 


(9) 


^i=l 


i=l 


Inequalities (|5]) and ([5]) are respectively known as the continuous-time Jensen inequality and the discrete-time 
Jensen inequality, which have been widely used in the literature for the stability analysis and stabilization 
of time-delay systems (see [5] and the references that have cited it). By using these two inequalities jointly 
we get the following corollary. 


Corollary 1 Let Ti > 0,i G I[l,fV], be N given nonnegative scalars. Assume that tVi : [— Ti,0] -G- R",i € 
I [I, TV], are such that the integrals in the following are well-defined, then 

/atn'^/wn, n 0 

(E^M *3 ( E^M - ^E'^*y wT’(s) Qwi (s) ds, (10) 

where Xi = Wi (s) ds, i G I [1, TV]. 
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We notice that all the N integrations ujf (s) QuJi (s) cls,i S I [1, N], on the right hand side of (|T0l) share 
the same weighting matrix Q, which is clearly very restrictive. To reduce the possible conservatism, we 
introduce the following multiple Jensen inequality. 


Lemma 4 Let Qi G R”^", i € I [1, A^], be N given positive definite matrices and Ti > 0,i G 1 [1, N], be N 
given scalars. Assume that the vector functions uji : [—ri,0] R",i G I [1,7V], are such that the integrals in 
the following are well-defined, then 

T 


N 


N 






(11) 




where Xi = uJi (s) ds, i € I [1, A^] and Q — Moreover, for a series of vectors € R^, i ^ I [1, A^] , 

there holds 


N 


N 


N 


N 


E^O E^O 




V2=l 


^2=1 


2=1 


Proof. Notice that, for any i ^ I [1, A^], by a Schur complement, there holds 

Taking integration on both sides of the above inequality gives 

1° (s) nQ-^u}, (s) ds f° ujJ (s) ds 


/_ u}^ (s) ds 


Qi 


> 0 , 


( 12 ) 


(13) 


(14) 


where 7 G I [1, TV] , which implies 


^ n ^ n 

E I-r, ‘^i (■S) nQf^UJi (s) ds J2 I-r, («) ds 
2=1 2=1 


N 


E f-r, (s) ds 


N 

EQ^ 
2 = 1 


> 0 . 


(15) 


By a Schur complement again, (HSl) is equivalent to CD- Finally, the inequality in (11211 can be proven in a 
similar way. ■ 

Now every integration (s) QiWi (s) ds,i G I [1,7V], on the right hand side of (ITT]) is weighted by 

an individual weighting matrix Qi,i G I[1,7V], which can introduce more decision variables that can be 
optimized to reduce the conservatism of the resulting conditions. The multiple Jensen inequalities m and 
(IT^ are clearly less conservative than the inequalities in (TTUI) and m since the later ones can be obtained 
immediately by setting Qi = Q~^,i G I [1, TV], in the former ones. 

By applying the Jensen inequality (1101) in Corollary [1] and choosing the following Lyapunov functional 


V (xt) = f {s) Px (s) ds 
J t—r 

N „0 

+ y] / (s + Tj) x’’" (7 + s) QjX (t + s) ds, 

.•_1 d— Ti 


(16) 


the following result was obtained in [Si- 


Lemma 5 176] / The IDS m is exponentially stable if there exist TV + 1 positive definite matrices P, Qi G 
R"^",7 G I[l, A^], such that the following coupled LMIs are satisfied 


N 


NnAj I -P + E I A, - Q, < 0, 7 G I [1, Af]. 

7=1 


(17) 


In the next two subsections, we will show how to use the multiple Jensen inequality (EH) to improve the 
above result. 
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3.2 The First Sufficient Stability Condition 

In this subsection we present a new sufficient condition for the exponential stability of the IDS ([T|) by applying 
the multiple Jensen inequality (El) and choosing a similar Lyapunov functional as (El- 

Theorem 1 Consider the IDS flj)- Then 


1. It is exponentially stable if there exist 2N positive definite matrices Si,Qi S G I[1,A^], such 

that the following nonlinear matrix inequalities are satisfied 


N / N \ 

'^Si < I Qi 


-1 


(18) 

(19) 


i=l 


v,i=l 


2. The nonlinear matrix inequalities in are solvable if and only if there exist 2N positive definite 

matrices Si, Qi G R"^", i G I [1, A^], and a matrix R G R"-^" such that the following LMIs are satisfied 


'y ' Qi + 'y ^ Si — + r) < 0 , 

i=l i=l 

-S, nAjR 

TiR^Ai -Qi 


< 0, f G I[l,iV]. 


( 20 ) 

( 21 ) 


Proof. Proof of Item 1. The inequality in (fT9|) implies that there exists a positive definite matrix P > 0 
such that P < eln for some sufficiently small number e > 0 and such that 


N 


N 




( 22 ) 


Consider the following Lyapunov functional 


V (xt) = f x"^ (s) Px (s) ds 

J t — T 

^ f s \ 

+ y] / 1 -h 1 I (t + s) SiX [t + s) ds, 

i=i “'-n / 


(23) 


which is in the form of (jl 6 |) that was used in El- The time-derivative of V (xt) satisfies 


N 


N 


V (xt) < x^ (t) I p+5^ j X (^) - 

\ 

C N N 

^ Qi I x{t)-'^yi 

N 0 

- E / (t + s) nAjQf^AiX {t + s)ds-^ y, 

i=i i=l 

^ 1 

= y] — / x'^ {t + s) (xiAjQf^Ai - Si) x{t + s) ds 

. , Xi J _T 


N 


< -7 


fit- 


‘ ds, 


(24) 


where Pi = ^ {t -\- s) SiX (t + s) ds, 7 > 0 is some constant and we have used the multiple Jensen 

inequality m and the nonlinear matrix inequalities (IT51) . Hence, it follows from Lemma [T] that the IDS (fi]) 
is exponentially stable. 
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Proof of Item 2. Let S — Q — J2^iQj- Then (fTOl) is equivalent to Q < S' which, by Lemma 

m is satisfied if and only if there exists an i? G such that 

R^QR + S - {R^ + R) < 0. (25) 

On the other hand, by the Schur complement, the inequalities in (1181) are satisfied if and only if 


-Si nAj 
'^iAi Qi 


<0, iei[i,iv]. 


(26) 


which, by a congruence transformation, are equivalent to 


-Si nAjR 

nR^Ai -R^QiR 


< 0, i G I[1,1V]. 


(27) 


It is clear that (1^ and (E7)) are respectively equivalent to (pH)) and (^11) by a substitution R^QiR Qi, i G 
I [1, N] (and thus R^QR Q). The proof is finished. ■ 


3.3 The Second Sufficient Stability Condition 

With the help of the multiple Jensen inequality (fill) , we further present in this subsection a new sufficient 
condition for the exponential stability of the IDS ([T|) with an alternative Lyapunov functional, which may 
possess some advantages over (Hi and (|23ll . 

Theorem 2 Consider the IDS (PP- Then 


1 . 


It is exponentially stable if there exist N positive definite matrices Qi G G I , such that 

the following nonlinear matrix inequality is satisfied 


N 

E 


r^AjQ-^Ai-(j2Q^ 


< 0 . 




(28) 


2. The nonlinear matrix inequality S28\} is solvable if and only if there exist N positive definite matrices 
Qi G R”^", i G I [1, iV] , such that the following LMI is satisfied 


N 

TlAl 


TlAl 

T 

Qi 

E 

2=1 

tnAn 

Q^ 

tnAn 

— 

Qn 


Proof. Proof of Item 1. Let Q = Then it follows from (l28t that there exist two sufficiently small 

numbers 6 > 0 and e > 0 such that 


N 

E + DSIn) - Q~^ < -eQ~^. 

2=1 


Let > 0, i e I [1, A^], be such that 


N 


N 


Y^r^Ar = q-^= 


and consider an associated nonnegative functional 

N 


u (xt) = E / 

i=l ■'t-Ti 


(30) 


(31) 


(32) 
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whose time-derivative is given by 


N 


N 


Vi {xt) = (t) I ^ j X (t) -^Vi 

/ AT \ N 

= (t) ( ) 2; (t) - ^ y. 




z=l 


N 


N 


< y] / T,x'^ {t + s) (i + s) ds - ^ yi. 


( 33 ) 


2=1 ' 


2=1 


where = x (t — Ti)^ RiX {t — ti) , and we have used the IDS H]) and the multiple Jensen inequality (ED- 
Choose another nonnegative functional 


^ rTi rt 


V 2 {xt) = '^ / x'^ {l){nAjQAA^ + SIn)x{l)dlds, 


/Q Jt—s 

whose time-derivative can be evaluated as 

N 

V 2 (xt) = (t) {TfAjQ~^A, + SrJri) x (t) 

2=1 

N .0 


/ u 

x^ (t + s) [nA^Q~^Ai + 5In) x(t + s) ds. 

-Ti 

2=1 ^ 


Hence, it follows from ([55|) and (I55]l that 

Vi {xt) + V 2 (xt) 


N 


N 


< (t) (xfAjQ^ ^A, + Snln) x(t) - '^yi - fj. 

2=1 2=1 
N 

< (1 - e) (t) Q~^x (t) - '^yi - fj, 

2=1 

N 

< (1 - e) y] (x'^ (t) RiX (t) -x(t- Ti)'^ RiX (t - - y 

2=1 

N „o 

= {l-e)Vi{xt)-S^ ||a;(t + s)f ds, 

A — l J-Ti 


N 


where y. = S f-r + ®)ll ds and we have used (1501) . Therefore 

i=l 


;(t ■ 


' ds. 


V (xt) = eVi (xt) + V 2 (xt) ^ J 

Finally, it is trivial to show that V (xt) satisfies ([5]). The conclusion then follows from Lemma [1] 
Proof of Rem 2. By a Schur complement, the LMI in (1001) is satisfied if and only if 


-Q 

tiQAJ ■ 

• xjifQAjf 

TiAiQ 

—Qi 


tnAnQ 


—Qn 


< 0 . 


By using a Schur complement again, the inequality (I38p can be equivalently transformed into 

N 


— Q + x^QAJQ i ^AiQ < 0 , 


(34) 


(35) 


(36) 


(37) 


(38) 


(39) 
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which is further equivalent to (E51) . The proof is finished. ■ 

At the end of this section, we will show that the stability condition (0^ has a very interesting relationship 
with the stability condition of the following IDS 


N 

X (t) ='^ A^x {t - Ti) , (40) 

i=l 

which was originally studied in [T] by using a Lyapunov functional approach, where G G I [1, A^], 

are given matrices and t^, i € I [1, iV], are given scalars and are such that 0 < n < • • • < tat. 


Lemma 6 0/ The IDS is exponentially stable independent of the delays Ti,i G I[1,A^], if there exist 
N positive definite matrix Qi G R"^", i G I [1, A^] , such that the following LMI is satisfied 


N 

■ Ai ' 


■ Ai ■ 

T 

Qi 

E 

Ajsf 


An 

— 

Qn 


Proof. By recognizing the characteristic polynomial of the IDS (HiH) (see Eq. (3.7) in [T]), it is not hard to 
see that it is exponentially stable if and only if the following IDS 

N 

x{t) = ^Ajx{t-Ti), (42) 

i=l 


is. Then by the results in Section 4 in [T], the IDS (1421) is exponentially stable independent of the delays 
Ti,i G I[l,Af] , if there exist N positive definite matrices Xi,i G I[l,Af] , such that the following LMI is 
satisfied (see inequality (4.6) in [T]) 


Hi -AiXiAj 

-A2X1AJ Ea 


—AiXiAjj 

—A2XiAjf 


> 0 , 


—Af^XiA^ —AatATiAJ 


Xn — AatATiA^ 


(43) 


where E^ = — AiXiAf — Xi+i,i S I [I, A — I]. As the above LMI implies Xi > Xi+i,i S I [1, A — 1], we 

can let Qat = Aat and Qi = Xi — >0, i G I [1, A — 1]. Then the LMI in (H51) can be exactly rewritten 

as (HD). 


It is very interesting to notice that the LMI (1291) and the LMI (|4T]) possess the very similar structures and 
the only difference is that the former LMI is delay dependent and the later one is not. This similarity may 
help us to understand the stability of these two classes of IDSs © and (HOI) 


4 A Spectral Radius Based Condition and A Comparison 

4.1 A Spectral Radius Based Stability Condition 

In this subsection, we will present a spectral radius based sufficient condition for the exponential stability of 
the IDS © based on Lemma [5l as indicated by the following theorem. 

Theorem 3 The following statements are eguivalent. 

A. There exist A + 1 positive definite matrices P, Qi G R""^", i G I [I, A], such that the coupled LMIs in 
are satisfied. 










B. There exist N positive definite matrices Qi S R"^”, z G I [1, A^], such that the following coupled LMIs 
are satisfied 

/ N \ 

NrfAj I ^ Q,1 >1, - Q, < 0, z G I [1, iV]. (44) 

C. The following condition is met, where p (A) denotes the spectral radius of a square matrix A : 

(45) 

D. There exists a positive definite matrix Q G such that the following LMI is satisfied 

N 

J2^rfAjQA,-Q<0. (46) 

i=l 


Proof. Proof of A-i^B. It is clear that the LMIs in (ITT)) are satisfied if we choose P = eln with £ being 
sufficiently small and the following coupled LMIs 

N \ 

A-Q,<0, i€l[l,N], (47) 

are satisfied. The converse is obvious. Finally, the LMIs (ITTl) are equivalent to (HU) by the substitution 
TiQi Qi, z G I [1, N]. 

Proof of B^C. Let 5”'^” = {Si,S 2 , ■■■ ,Sn) where Si G ^ {S : S = > 0}. Then 



-^(Q) 


N 

1=1 


N \ 

■■ , NrfjAjf QjAjn I , 


(48) 


where Q = {Qi, ■ ■ ■ ,Qn) G is a linear positive operator (see Definition 1 in [H]). Consequently, the 

inequalities in (HU) are satisfied if and only if there exists a Q G 5"^" such that 

^{Q)-Q<0, (49) 


where V < 0 means —V G Then by Lemma I in the inequality in (H^ has a solution Q G 

if and only if p (J£') < 1. However, similar to (lO)-(II) in [IT], we can show that p (.if) = p = p {NA) 

with 


rfAj ® Aj rfAj ® Aj 

rfA^ (g) A^ ® A^ 



rfAj ® Aj 

0 A^ 

0 Ajj 


Hence the LMIs in (H^ are solvable if and only if 


p{A) < 


N 


Now notice that we can write A = BC where 

rfAj 0 Aj 


B = 


t%A^ 0 


,C=[ I 




(50) 


(51) 


On the other hand, for any two matrices X and Y with appropriate dimensions, we have p {XY) 
Hence 


p{A)=p{BC) 


piCB)=p 



p{YX). 


(52) 
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by which the inequality in (15(11) is exactly the one in (H51) . 

Proof of C^D. The proof is quite similar to the proof of by utilizing the linear positive operator 

N 

^{Q) = Y.NrfAjQA, (53) 

i=l 

where Q S S"^". The proof is finished. ■ 

Item C of Theorem [3] implies an interesting spectral radius based sufficient condition for testing the stability 
of the IDS ([1]), as highlighted in the following corollary. 

Corollary 2 The IDS m is exponentially stable if the spectral radius condition is satisfied. Particu¬ 
larly, if N = 1, the IDS m is exponentially stable if p{Ai) < -A. 

Remark 1 If N = 1, by Proposition 2 in m, the IDS m is exponentially stable if ||Ai|| < which is 
more conservative than p{Ai) < ^ since p{Ai) < ||Ai|| for any matrix Ai. 

Corollary 3 The IDS m is exponentially stable if there exist N scalars at G (0,1),* G I[l,iV] such that 
= 1 and 

P < 1. (54) 


Proof. According to the proof of Theorem [3l (|54|) is satisfied if and only if there exists a Q > 0 such that 

N 2 

Y^AJQA,-Q<0. (55) 

Hence the inequality (1^51) is satisfied with Qi = UiQ. The result then follows from Theorem [5] ■ 

Clearly, the spectral condition (15^ reduces to (l45ll if we set Oi = * G I [1, A^]. 

4.2 A Comparison of Different Sufficient Conditions 

With the help of Theorem [31 we are able to make a comparison among these different stability conditions in 
Lemma [51 Theorem [I] and Theorem [31 


Proposition 1 The following statements are true. 

1. If the set of LMIs m ( or M) are solvable, then the set of LMIs in are also solvable, 

namely. Theorem [7] is always less conservative than Lemma O 

2. The set of LMIs in \20\) -i21 \) are solvable if and only if the LMI in i2ffl) is solvable. Hence Theorem 
H is equivalent to Theorem[J\ and both of them are thus less conservative than Lemma O 


Proof. Proof of Item 1. By Theorem [31 the set of LMIs (1171) are solvable if and only if the set of LMIs (l44l) 
are solvable, which implies that there exists a sufficiently small number e > 0 such that £/„ < Qi,i € I[l, N], 
and 

/ N \ 


NrfAj 


i=i 


Qj I A ~ Qi + sin <0, * G I [1, N] . 


(56) 
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Now for every i € I [1, A^] we choose 


N 


p-^ = NY,Q^ = NQ, R^ = Q^- eln > 0, 


(57) 


which implies 


N 


N 




N 


Then it follows from ((56ll that, for all i S I [1, A^], 

T^AJP-^A, -R, = Nr^AjQAi - Q, + e/„ < 0, 

and it follows from (i58l) that 


N 


N 


N 


-1 


(58) 


(59) 


(60) 


Y^R,=J2{Q,-£l^) = Q-NeIn<Q= ■ 

i=l 1=1 \i=l / 

Clearly, (l59l) and (1601) are respectively in the form of (fT8l) - (|T^ . 

Proof of Rem 2. Assume that (|18|) " (fTOl) are feasible. Summing the N nonlinear matrix inequalities in dT51) 
on both sides gives 

N N 

Yr?AjQ-^A,-YS^<0, (61) 

which, by using the nonlinear matrix inequality f|19D . implies 

-1 


N 


N 


'A < 


(62) 


which is just in the form of f|28p and is further equivalent to 
Now assume that (l28l) is feasible. Denote 

N \ N 






2=1 


and let 


S, = T^AjQ-^A, + n>0, iGl[l,iV]. 
It follows that (HU) is satisfied. Now, by (^51) , we have 


N 


N 


N 


N 


E^* = o E^?d + < E^?» 


2 = 1 


^ 2=1 


2=1 


^ 2=1 


(63) 

(64) 


which implies that (1191) is satisfied. The proof is finished. ■ 

This proposition demonstrates in theory that the multiple Jensen inequality m used in the proofs of 
Theorems [1] and [2] can indeed reduce the conservatism of the resulting stability conditions. 

Remark 2 Though Theorems [7] and are equivalent hy Proposition [U Theorem obtained by the novel 
Lyapunov functionals and |g/(| ) possesses an advantage over TheoremUi To see this, we notice that the 
total row size (denoted by and the total number of scalar decision variables (denoted by'^) in the LMIs 
of TheoremsUl and\^ are, respectively, given by 


‘I’Th.i = n + 2n, 'I'Th.i = n (n + 1) TV + 


4>Th.2 = riN, 'I'Th.2 = 


n(n+l) 


N. 


(65) 


R is well known that the computational complexity of an LMI is bounded by where p, is a constant 

(see Ii22j/). Hence the computation complexity of the LMIs in Theorem\^is significantly lower than that in 
T/ieorem [7J 
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Combining Lemma IHl Proposition [T] and Theorem [3] gives the following corollary. 


Corollary 4 The IDS |^d| j is exponentially stable independent of the delays Ti,i Gl [1, N], if p ® j < 

N- 


Table 1: The maximal allowable T 2 by using different approaches 



LMI dlB 

LMIs (HZl) 

LMI dMl) 

Condition (I45|l 

n 

= 0.4 

0.0317 

infeasible 

infeasible 

infeasible 

n 

= 0.3 

0.1146 

0.0474 

0.0474 

0.0474 

n 

= 0.2 

0.2418 

0.1527 

0.1527 

0.1527 

Tl 

= 0.1 

0.4882 

0.3414 

0.3414 

0.3414 


5 A Numerical Example 


We consider an IDS in the form of © with two delays, namely, 


x{t) = Ai [ x{t + s) ds + A 2 

J -Ti 


X ft + s) ds, 


( 66 ) 


where 


Ai 


-4 1 

-13 2 


0 - 
1 0 


(67) 


If A 2 = 0, this system has been studied in m- It is shown there that the LMIs in Lemma [5] are feasible if 
and only if 0 < ri < 0.4473 = r*. On the other hand, direct computation gives p{TfAi) = 0.9999, which 
clearly validates Corollary [21 


For a fixed ri, the maximal value of T 2 such that the LMI in (l29ll . the LMIs in (HB, the LMI in (l46ll . and 
the spectral condition (HSl) are feasible can be respectively computed by a bisection method. The results 
are recorded in Table 1. From this table we can observe that Theorem [2] is always less conservative than 
Lemma [5] established in m, which indicates that our approach based on the multiple Jensen inequality 
can considerably reduce the conservatism in the stability analysis of this class of IDSs. Moreover, the LMIs 
CB, the LMI (j46ll and the spectral condition (1451) lead to the same result, which validates Theorem |3l We 
mention that the results obtained by Theorem |T] are the same as that by Theorem |21 However, from the 
computational point of view. Theorem [3] is recommended to use in practice as it only involves one constraint 
and a single decision variable. 


To illustrate Corollary|3l let ti = 0.4 and T 2 = 0.02. From Table I we can see that (1451) is not satisfied. Since 
||Hi|| ^ ||^2||, we may let the weighting factor ^ of {tiAi) 0 {tiAi) be small enough so that the spectral 
radius of the resulting matrix {DAi) 0 (TiAi) is less than one. Indeed, if we choose ai = 0.9 and 

02 = 0.1 we can compute p(Ei=i^ ijiAi) 0 (xiAi)) = 0.9783, which implies the asymptotic stability of IDS 
(EB in this case. 


We finally mention that the conclusions obtained in the above have been approved by many other randomly 
chosen numerical examples that are not included here to save spaces. 


6 Conclusion 

This note has studied the stability analysis of a class of integral delay systems (IDSs) with multiple delays, 
which have wide applications in the stability analysis of time-delay systems, especially for neutral time- 
delay systems. By generalizing the well-known Jensen inequality to the case with multiple terms through 
introducing multiple weighting matrices, two Lyapunov functional based approaches have been established 
to yield set of sufficient stability conditions. Moreover, it is shown by the positive operator theory that 
the obtained new conditions are always less conservative than the existing ones and a spectral radius based 
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sufficient condition is obtained simultaneously. A numerical example has demonstrated the effectiveness of 
the established approaches. 
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